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Abstract 

We present Monte Carlo simulations of hysteresis loops of a model of a magnetic nanoparticle with a ferromagnetic 
core and an antiferromegnetic shell with varying values of the core/shell interface exchange coupling which aim to 
clarify the microscopic origin of exchange bias observed experimentally. We have found loops shifts in the field direction 
as well as displacements along the magnetization axis that increase in magnitude when increasing the interfacial 
exchange coupling. Ovelap functions computed from the spin configurations along the loops have been computed to 
explain the origin and magnitude of these features microscopically. 
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1. Introduction 

Compound magnetic structures formed by a ferro- 
magnet (FM) in close contact with an antiferromag- 
net (AFM) show interesting proximity effects which 
are enhanced when going to the nanoscale range [1] . 
Among them, one of the most studied phenomenon 
is the displacement of the hysteresis loops along 
the field axis observed in layered FM/AFM systems 
and core/shell nanoparticles when measured after 
cooling in a magnetic field, an effect that has been 
termed exchange bias (EB) after its first observa- 
tion by Meiklejohn and Bean in Co particles [2]. 
In spite of the intensive research developed in the 
last decades, there are still many open issues to be 
understood and no clear-cut model has been able 
to account for all the observed phenomenology [3]. 
Several recent experiments based on spectrocopic 
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techniques [4] have focused on the observation of 
the magnetization reversal at the FM/AFM inter- 
face in different layered systems since it is widely 
accepted that the magnitude of the exchange bias 
field is related to the existence of uncompensated 
moments at the interface. However, in nanoparticles 
with FM core/ AFM shell structure, details about 
the core/shell interface cannot be accessed by these 
techniques and only indirect information about the 
magnetic order at the inerface can be gained by mag- 
netization measurements. 

In this study, we present the results of Monte 
Carlo simulations of hysteresis loops of an indi- 
vidual core/shell nanoparticles, which, apart from 
the usual shift in the field axis direction, display 
asymmetry between the positive and negative field 
branches and also shifts in the vertical direction 
as observed in some experimental systems [5]. By 
inspection of the spin configurations during field 
cycling, we will show how this phenomenolgy is 
related to the peculiar interplay between exchange 
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coupling at the core/shell interface and the mag- 
netic order of its uncompensated spins induced after 
the field cooling process. 

2. Model and Results 

The details of model for the core/shell particle on 
which simulations are based were already presented 
in our previous studies [6], where the details about 
the simulation protocol used for the Monte Carlo 
method can also be found. The simulations are based 
on the following Hamiltonian: 

— = — Jij Si ■ Sj — ki (S*) 2 — h ■ Si , 

(i,3) ' « 

where Si are classical Heisenberg spins of unit mag- 
nitude placed at the nodes of a sc lattice. In the first 
term, the value of the n.n. exchange constants Jij 
depends on the spins belonging to different particle 
regions. At the core, Jy is FM and has been fixed 
to Jc = +10 K. AF spins at the shell have Jsh = 
— 0.5 Jq. The exchange coupling at the interface [de- 
fined as those spinson the core (shell) with at least 
one neighbor on the shell(core)] is Ji n t- 

Hysteresis loops obtained for a particle with to- 
tal radius R = 12a, surface shell thickness i?sh = 
3a and Jsh = — 0.5Jc, obtained after cooling from 
above the Neel temperature of the AF shell in a field 
hpc — 4 are shown in Fig. 1 (circles)for three values 
of Jint- As can be clearly seen in the upper panels 
of the figure, the loops display an increasing shift 
towards the negative field direction with increasing 
values of Ji n t, demonstrating that the observance of 
EB is related to the AF coupling of the shell spins 
the core spins at the interface. 

A detailed inspection of the microscopic config- 
urations attained after the FC process [6] revealed 
that, inspite of the AF character of the coupling be- 
tween the shell spins and also between the interfacial 
spins, a net magnetization component along the field 
direction is induced at the shell interface due to the 
geometrical symmetry breaking and the alignment 
of groups of spins into the field direction. This net 
magnetic moment generates local fields on the core 
spins that point into the same direction as the ex- 
ternal field, causing the shift of the hysteresis loops. 
This is further supported by the observation that 
the hysteresis loops are shifted by the same amount 
but towards the positive field axis when obtained af- 
ter cooling in a field applied in a direction negative 
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Fig. 1. (Color online) Upper panels display the hysteresis 
loops for a particle with radius R = 12 a obtained after FC 
down to T = 0.05 in a field h-pc — 4 (circles) and ftpc = —4 
(dashed lines) for three values of the exchange coupling con- 
stant Ji n t at the core/shell interface. Lower panels show the 
average magnetization projection of the core spins along the 
field axis (squares) and the hysteresis loops for the com- 
ponent of the magnetization transverse to the field direction 
Mtr (circles). 

with respect to the measuring field (see for example 
the dashed lines in Fig 1. for Jint/Jc = —0.5,-1). 

In addition, a clear asymmetry between the upper 
and lower loop branches developes with increas- 
ing values of the interface coupling. This feature 
can be more clearly seen by looking at the aver- 
age absolute value of the magnetization projection 
along the field axis through the reversal process, 
= \Si • z\ displayed in the lower panels of 
Fig. 1 for the core spins. The origin of the loop 
asymmetry can be traced by monitoring the values 
the so-called overlap q(h) and link overlap qh{h) 
functions along the hysteresis loops, that are a gen- 
eralization of similar quantities commonly used in 
the spin-glass literature [7] and that are defined 
as q(h) = jfJ2 i=1 Si(h F c) ■ S t (h) and q L (h) = 
lv7 Si(h FC ) ■ Sj(h FC ) S*i(h) ■ Sj(h), where in 
qh{h) the summation is over nearest neighbors and 
Ni is a normalization factor that counts the num- 
ber of bonds. An example of the field dependence 
of these overlaps computed only for the interfacial 
spins is shown in Fig. 2, where we have separated 
the contribution of the shell and core spins. A de- 
parture of qh from unity is known to be proportional 
to the surface of reversed domains formed at field h 
and, therefore, q^ is sensible to the existence of do- 
main walls. The sharp decrease of q^ for core spins 
and the symmetry of the peak around the negative 
coercive field indicates uniform reversal. However, 



2 



- 

cr 



(a)< 



-4-3-2-101234 
h(K) 




-4-3-2-10123 4 

h(K) 



Fig. 2. (Color online) Field dependence of the link overlap 5l 
(a) and overlap q (b) functions for the interfacial spins at the 
shell (circles) and at the core (squares) for Ji n t = — 0.5Jc- 

the progressive reduction of along the ascend- 
ing branch and the asymmetry of the peak around 
the positive coercive field indicates the formation 
of domain walls at the particle core that sweep the 
particle during reversal. The function qih) measures 
differences of the spin configuration at field h with 
respect to the one attained after FC. Therefore, 
the decrease of q for the interface shell spins when 
reducing the magnetic field indicates the existence 
of a fraction of shell spins that reverse dragged by 
core spins, while the constancy of q in the ascending 
branch reveals the existence of spins pinned during 
core reversal. 

Clearly correlated to the observation of EB and 
the loop asymmetry, the loops experience a shift 
along the vertical (M z ) axis which increases with 
Jint, as reflected in Fig. 1 by the difference of the 
M z values in the high field region of the two loop 
branches or at the remanence points. The field de- 
pendence of magnetization transverse to the field 
direction M tr (circles in the lower panels of Fig. 
1), indicates that M tr has higher values for the de- 
scending loop branch that in the ascending branch, 
and that M tr increases with J\ n t. Snapshots of the 
spin configurations at the remanence points shown 
in Fig. 3, show the existence of a higher amount of 
core spins with transverse orientation near the in- 
terface at the lower branch (Fig. 3b, d) than at the 
upper branch. This is in agreementwith the results 
of some experiments of oxidized particles where this 
vertical shift was also observed [5]. Our simulation 
results above, indicate that the microscopic origin 
of the vertical shift is the due to the different rever- 
sal mechanisms on the two loop branches due to the 
existence of uncompensated pinned moments at the 
core/shell interface that faciltate the nucleation of 
non-uniform magnetic structures during the ascend- 
ing field branch of the loops. 

We acknowledge H. Katzgraber for helpful 
comments on overlap functions and CESCA and 
CEPBA under coordination of C 4 for computer fa- 
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Fig. 3. (Color online) Snapshots of the remanent spin con- 
figurations of the upper (a,c panels) and lower (b,d panels) 
branches of the hystersis loops showing midplane cross sec- 
tions of the particle parallel (a,b) and perpendicular to the 
z axis (c, d) for the case Ji n t = — Jc shown in Fig. la. Dark 
(light) blue cones represent core (core/interface) spins while 
green (yellow) ones are for spins at the shell (shell/interface). 

cilities. This work has been supported by the Span- 
ish MEyC through the MAT2003-01124 project 
and the Generalitat de Catalunya through the 
2005SGR00969 CIRIT project. 



References 

[1] J. Nogues et al., Phys. Rep. 422 (2006) 65. 

[2] W. H. Meiklejohn and C. P. Bean, Phys. Rev. 102 (1956) 
1413; Phys. Rev. 105 (1957) 904. 

[3] M. Kiwi, J. Magn. Magn. Mat. 234 (2001) 584; R. L. 
Stamps, J. Phys. D: Appl. Phys. 33 (2000) R247. 

[4] H. Ohldag et al, Phys. Rev. Lett. 87 (2001) 247201 
H. Ohldag et al, Phys. Rev. Lett. 91 (2003) 017203 
P. Kappenberger et al., ibid. 91 (2003) 267202 
P. Blomqvist et al., ibid. 94 (2005) 10720.3; S. Roy et al, 
ibid. 95 (2005) 047201; V. K. Valev et al., ibid. 96 (2006) 
067206; E. Arenholz et al, Appl. Phys. Lett. 88 (2006) 
072503. 

[5] R. K. Zheng et al, Phys. Rev. B 69 (2004) 214431; 
L. Del Bianco et al, ibid. 70 (2004) 052401; J. B. Tracy 
et al, ibid. 72 (2005) 064404; E. Passamani et al, J. 
Magn. Magn. Mat. 299 (2006) 11. 

O. Iglesias et al, Phys. Rev. B 72 (2005) 212401; 
O. Iglesias and A. Labarta, Physica B 372 (2005) 247. 



[6] 
[7] 



H. G. Katzgraber et al, Phys. Rev. Lett. 89 (2002) 
257202; Phys. Rev. B 63 (2001) 184422. 



3 



